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Abstract 

In this work, the charm quark mass is obtained from a QCD sum rule 
analysis of the charmonium system. In our investigation we include 
results from nonrelativistic QCD at next-to-next-to-leading order. Us- 
ing the pole mass scheme, we obtain a value of M c = 1.70 ± 0.13 GeV 
for the charm pole mass. The introduction of a potential-subtracted 
mass leads to an improved scale dependence. The running MS-mass 
is then determined to be m c {m c ) = 1.23 db 0.09 GeV. 
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1 Introduction 



Quantum Chromo Dynamics (QCD), describing the strong interactions, is 
one of the key components of the Standard Model. The determination of its 
parameters remains an essential task within modern particle physics. Much 
effort has therefore been put into the study of quark masses. However, in 
most cases where the system is sensitive to mass effects, also confinement 
plays an important role. Apart from the top mass, the quark masses can thus 
only be calculated by nonperturbative methods like QCD sum rules pHl, 
lattice QCD ||,||] or chiral perturbation theory f|[7|]- 

In the past, QCD moment sum rule analyses have been successfully ap- 
plied for extracting the charm and bottom quark masses from experimental 
data on the charmonium and bottonium systems respectively. The funda- 
mental quantity in these investigations is the vacuum polarisation function 
ri(g 2 ): 

IV(g 2 ) = i Jd'x e 4 *" (T{jp(x)j u (0)}) = (q^ - g, v q 2 ) I% 2 ) , (1) 

where in the charm case the vector current is represented by = (c7 M c)(x). 
Via the optical theorem, the experimental cross section a(e + e~ —>■ cc) is re- 
lated to the imaginary part of U.(s): 

„, . 1 a(e + e~ — > cc) , 

Using a dispersion relation, we can express the moments A4 n of II(s) by 
integrals over the spectral density R(s): 



12vr 2 / 2 d\ n , , 



^rr^- (3) 

S=0 J S m in * 



For convenience, we have defined the moments as dimensionless quantities. 
In addition, it will prove useful to express the moments by integrals over the 



velocity v = Jl — Am? j s: 



M n = 2 dvv(l - v 2 ) n - L R(v) . (4) 







The moments can either be calculated theoretically, including perturba- 
tion theory, Coulomb resummation and nonperturbative contributions, or 
can be obtained from experiment. In this way, we can relate the charm 
quark mass to the hadronic properties of the charmonium system. 
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During the last years, several analyses for the bottonium system have 
been performed [p Hl5|1 . In these investigations, it turned out that the largest 
theoretical contributions arise from the threshold expansion in the framework 
of nonrelativistic QCD (NRQCD) which is known up to next-to-next-to- 
leading order (NNLO) |12], |14|,|16|1 . The correlator can be expressed in terms 
of a Greens function, which shows a continuous spectrum above threshold, 
but also contains poles below threshold. 

Yet these new theoretical developments have not been applied to the 
charmonium system. Due to the low scales involved, the analysis becomes 
more delicate here. We have therefore introduced a new method of analysis 
which allows for investigation of the influence of the individual contributions 
on the error. In particular, we have taken special care in obtaining the pole 
contributions and the reconstruction of the spectral density above threshold. 

So far, we have not specified the mass definition to be used in eq. ([IJ). 
A natural choice for this mass is the pole mass M. In the first part of our 
numerical analysis, we will use the pole mass scheme to extract the charm 
pole mass. However, as the pole mass suffers from renormalon ambiguities 



171 , it can om y be determined up to corrections of order Aqcd- In the second 



part of our analysis we shall therefore use the recently introduced potential- 
subtracted mass nips (HH- From this mass definition we can obtain the 



MS-mass more accurately than from the pole mass scheme. 

In the next section, we shall present the contributions from the threshold 
expansion in the framework of NRQCD. The perturbative expansion which 
is needed for the reconstruction of the spectral density will be discussed in 
section 3. Afterwards, we will shortly describe the nonperturbative contribu- 
tions and the phenomenological spectral function. In the numerical analysis, 
we shall first explain the method of analysis. We will then obtain the pole 
mass and the MS-mass from an analysis of the pole mass and the PS-mass 
scheme respectively. The origin of different contributions to the error will be 
carefully investigated. We shall conclude with a summary and an outlook. 



2 Coulomb resummation 

Close to threshold, it turns out that the relevant expansion parameter is a s /v 
rather than a s since the velocity becomes of order a s . These terms can be 
resummed in the framework of NRQCD and give the largest contributions 
to the theoretical moments. 



The correlator is then expressed in terms of a Greens function [12, 19 



n(5) = 2^ [ C »M G W + m-2 G c(k)) , (5) 
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where k = y M 2 — s/4 and M represents the pole mass. The constant Ch(o> 
is a perturbative coefficient which is needed for the matching between the full 
and the nonrelativistic theory. It naturally depends on the hard scale. The 
Greens function is analytically known up to NNLO jn| and sums up terms of 
order a n s jv n ^ k for n > and k = 1,2, 3. It is crucial for the analysis that the 
result depends on three scales. While the hard scale fjthard is responsible for 
the hard perturbative processes, the soft scale \i so ft governs the expansion of 
the Greens function. Furthermore, the factorisation scale fif ac separates the 
contributions of large and small momenta and plays the role of an infrared 
cutoff. 

The Greens function contains two parts: the continuum above and the 
poles below threshold respectively. We are interested in both contributions 
separately: first, the individual corrections can be analysed and their error 
estimated. Second, in our numerical analysis we will reconstruct the spectral 
density above threshold and we thus need the corresponding spectral density 
at low velocities. In principle, the expressions for the energies and decay 
widths of the poles have been calculated. However, in the actual case of 
the charm quark the expansion does not converge well. We will therefore 
choose a different method of evaluation [P0| , |2Tfl . Since the expansion for the 
Greens function at NNLO is known analytically [0, we can evaluate their 



contribution to the moments numerically by performing the derivatives. On 
the other hand, by using a dispersion relation, we can derive the continuum 
from the imaginary part of the correlator. From the difference we can then 
obtain the pole contributions: 



M Poles = ( AM 2~ 




2 x" r°° , ImII(s) 



127T UM'\ ds—^L. (6) 



4M 2 S n+1 



Since we will not evaluate the poles near threshold, but rather calculate 
their contributions to the moments in a region where perturbation theory is 
expected to be valid, the convergence of the pole contributions is improved. 
Nevertheless, the poles will give the largest contribution to the theoretical 
moments and thus the dependence on the scales will remain relatively strong. 
In the numerical analysis we will give a detailed account on the size and 
behaviour of these contributions. 

The large corrections are partly due to the definition of the pole mass. 
These contributions can be reduced by using an intermediate mass definition. 
In this analysis we will use the so-called potential-subtracted (PS) mass fT8|j 
where the potential below a separation scale fi sep is subtracted: 

1 f cPq 

mps(fisep) = M - 5m(fi sep ) , 5m(fi sep ) = -- J -^—^ V(q) . (7) 

|q|<Msep 
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Since the QCD potential contains the same renormalon ambiguities as the 
pole mass, in the difference they are cancelled. This mass definition thus 
leads to an improved scale dependence and a more precise determination of 
the MS-mass. 



3 Perturbative expansion 

The perturbative spectral function R(s) can be expanded in powers of the 
strong coupling constant a s , 

RP\s) = R®(s) + ^ RM(s) + ^ R {2) (s) + .... (8) 

From this expression the corresponding moments M. 1 ^ can be calculated via 
the integral of eq. (f|). The first two terms are known analytically and can 
for example be found in ref. [§]. I^ 2 \s) is still not fully known analytically 
However, the method of Pade-approximants has been exploited to calculate 
n^ 2 ) numerically, using available results at high energies, analytical results 
for the first eight moments and the known threshold behaviour p2| , |23[] . This 
information is sufficient to obtain a numerical approximation of Yi^ 2 \s) in 
the full energy range. 

The numerical stability of the results can be checked in different ways. By 
choosing different Pade-approximants or by selecting a smaller set of input 
data the results for the moments remain almost unchanged. Furthermore, 
some contributions to the spectral density like those from internal quark 
loops are known analytically and are in very good agreement with the nu- 
merical spectral density. This provides strong support that the numerically 
obtained spectral density gives a very good approximation to the exact spec- 
tral density. 



4 Condensate contributions 



The nonperturbative effects on the vacuum correlator are parametrised by 
the condensates. The leading correction is the gluon condensate contribu- 
tion which is known up to next-to-lcading order |24|. Furthermore, the di- 



mension 6 and 8 contributions have been calculated and will be included in 
our analysis [p5|-27|. In the analysis below, we have employed a value of 
(a s /7rFF) = 0.024 ± 0.012 GeV 4 for the gluon condensate. 

It will turn out that the condensate contributions are suppressed when 
compared to former charmonium sum rule analyses [^,^] and only have 
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little influence on the mass. Besides an increase of the theoretical moments 
from the Coulomb contributions we will restrict the moments to n < 7 where 
the nonperturbative contributions are relatively small. Since we obtain a 
larger pole mass than the former analyses, the condensates, starting with a 
power of 1/M 4 , are suppressed further. 



5 Phenomenological spectral function 

Experimentally, the lowest lying six ^-resonances have been observed. Since 
the widths of the poles are very small compared to the masses, the narrow- 
width approximation provides an excellent description of these states. To 
model the contributions of the hadronic continuum, we use the assumption 
of quark-hadron-duality and integrate over the perturbative spectral density: 

To estimate the continuum contribution we will use a threshold in the range of 
3.6 GeV < yfso < 4.2 GeV with a central value of = 3.8 GeV. The lower 
value corresponds to using a perturbative continuum around the ip(2S) but 
nevertheless including the higher resonance contributions. This is certainly 
over-counting the continuum contribution. The higher value would assume 
no additional continuum contribution until the ip{5S) resonance. However, 
the most dominant phenomenological contributions come from the first two 
^-resonances resulting in a small influence of the continuum even for low 
values of n. 



6 Numerical analysis 
6.1 Method of analysis 

Besides the contributions from the poles of the Greens function and the con- 
densates, the theoretical part of the correlator contains the spectral density 
above threshold. Now we will discuss the different parts of the spectral den- 
sity. A more detailed description of our procedure will be presented in a 
forthcoming publication . 

For high velocities the spectral density is well described by the perturba- 
tive expansion where we have fixed the hard scale to the pole mass. However, 
as one approaches smaller values of v, the perturbative expansion breaks 
down. The resummed spectral density, on the other hand, gives a good de- 
scription for low values of v, but becomes unreliable for high velocities. For 
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Figure 1: Thick solid line: reconstructed spectral density; Thin solid line: 
perturbative spectral density; dotted line: perturbation theory at NLO; 
dashed line: resummed spectral density. 



these reasons, we will introduce a separation velocity v sep ~ 0.3. Above 
v sep we will use the perturbative spectral density. Below v sep , we take the 
resummed spectral density adding the terms which are included in perturba- 
tion theory but not in resummation. 

In fig. [I] we have displayed the different contributions. The dotted line 
represents the perturbative expansion at NLO. The thin solid line also in- 
cludes the NNLO. Whereas for high velocities the perturbative expansion is 
well convergent, the importance of the higher corrections increases for smaller 
v. The dashed line is the resummed spectral density and the thick solid line 
the reconstructed spectral density. For the charmonium system, there exists 
a range of intermediate values of v where neither the perturbative expansion 
nor the resummation can be trusted. Indeed, it can be clearly seen that the 
reconstructed spectral density shows a gap at the separation velocity. To 
estimate the error we have varied v sep between 0.2 and 0.4. The analysis 
shows that though the introduction of the separation velocity stabilises the 
sum rules, the variation only has a minor influence on the mass. 



6 



Mr 



1.9 - 




n 



Figure 2: Thick solid line: central pole mass; thin solid lines: M c for 



Vsoft — 1-05 and 1.5 GeV; dashed lines: M c for fif a 
dotted lines: M c for fihard = 1-4 and 2.5 GeV. 



1.2 and 1.7 GeV; 



6.2 Pole mass scheme 

Since the threshold region becomes more important for larger values of n, as 
can already be seen from the factor (1 — v 2 ) n ~ l of eq. (Q), also the NLO and 
NNLO corrections grow for larger n. We will therefore restrict our analysis to 
moments with n < 7. As our method of analysis needs moments of n > 3 to 
reconstruct the spectral density, we will use 3 < n < 7. Since the continuum 
part of the phenomenological spectral density is relatively small even for low 
values of n, it does not set a lower bound on n. As central values for our 
scales we have selected: 

fx soft = 1.2 GeV , fx fac = 1.45 GeV, /w = 1.7 GeV . (10) 

We have set the hard scale equal to the central value for the pole mass. For 
the soft scale, we would have liked to choose a value around ^ so ft ~ 0.8 — 0.9 
GeV, but the NNLO corrections get out of control for these scales. Thus, 
we take a central value of fi S oft = 1-2 GeV and vary it in an interval of 
1.05 GeV < fi so ft < 1.5 GeV. The factorisation scale separates the different 
regions and should lie between the two other scales. 

In table [1] we have collected the individual moments for different values 
of n. A4^ oles are the theoretical poles of the Greens function, j\z[^ RResum the 
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n 


3 


4 


5 


6 


7 


KA Poles 


3.84 


4.44 


5.11 


5.87 


6.73 


J l n 


1.92 


1.55 


1.32 


1.15 


1.03 


KA N RResum 


0.28 


0.27 


0.26 


0.24 


0.23 


KA Counter 
J l n 


0.69 


0.66 


0.64 


0.62 


0.60 


KA C ontinuum 
J l n 


0.75 


0.47 


0.31 


0.21 


0.15 


kaC ondensates 
J l n 


-0.03 


-0.05 


-0.06 


-0.08 


-0.09 



Table 1: Moments for different n with the parameters [i so ft — 1-2 GeV, 
fifac = 1.45 GeV, fi hard = 1.7 GeV, ^/s^ = 3.8 GeV and v sep = 0.3. 



Hsoft 
J^Poles 


1.05 
9.50 


1.1 
7.52 


1.2 
5.14 


1.3 
3.80 


1.4 
2.97 


1.5 
2.41 


Ufac 
J^Poles 


1.2 
6.79 


1.3 
6.12 


1.4 
5.46 


1.5 
4.82 


1.6 
4.18 


1.7 

3.55 


fihard 
J^Poles 


1.4 
4.24 


1.5 
4.62 


1.6 
4.93 


1.8 
5.45 


2.0 
5.85 


2.5 
6.56 



Table 2: Pole moments with n = 5 for different fi so ft with fif ac = 1.45 GeV 
and \ihard = 1.7 GeV, for different fif ac with (j, so ft = 1-2 GeV and [ihard = 1.7 
GeV and for different fihard with fj, so ft = 1-2 GeV and fj,f ac = 1.45 GeV. 



moments from the nonrelativistic spectral density and j\y[^ ounter contains the 
double counted terms from the perturbative expansion below v sep . In this 
scheme the theoretical moments are dominated by the pole contributions. 
Though the size of these terms grow with larger n, the mass remains relatively 
stable since the contributions of the moments to the mass is suppressed by 
a power of l/2n. 

The thick solid line in fig. |2| shows the central value for the pole mass 
M c = 1.70 GeV with the parameters from eq. fllUp. The error is dominated 
by the variation of the scales. We obtain: 

1.05 GeV < n soft < 1.5 GeV : AM C = 100 MeV , 

1.2 GeV < n fac < 1.7 GeV : AM C = 50 MeV , 

1.4 GeV < fihard < 2.5 GeV : AM C = 40 MeV . (11) 

In particular, in table |2| we have listed the influence of the scales on the 
pole contributions. We have varied the scales within physically reasonable 
ranges. Since the convergence of the nonrelativistic expansion is not very 
good for the charmonium system, the scales cannot be chosen arbitrarily 
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t^soft 




1.05 


1.1 


1.2 


1.3 


1.4 


1.5 




LO 


3.03 


2.57 


1.95 


1.56 


1.29 


1.10 


J^Poles 


NLO 


5.60 


4.75 


3.61 


2.89 


2.40 


2.04 




NNLO 


9.50 


7.52 


5.14 


3.80 


2.97 


2.41 




LO 


0.64 


0.61 


0.57 


0.54 


0.51 


0.49 


Ji^ N RResum 


NLO 


0.24 


0.27 


0.29 


0.30 


0.31 


0.31 




NNLO 


0.31 


0.28 


0.23 


0.20 


0.18 


0.17 



Table 3: Size of the moments from the poles and the resummed spectral 
density at LO, NLO and NNLO for different values of \i ao ft- 



Source 


AM C 


Variation of fM so ft 
Variation of /i/ ac 
Variation of fihard 
Threshold So 
Experimental error 
Variation of v sep 
Variation of Aqcd 


100 MeV 
50 MeV 
40 MeV 
10 MeV 
15 MeV 
10 MeV 
50 MeV 


Total error 


130 MeV 



Table 4: Single contributions to the error of M c . 



far away from their central values. Though the analysis is stable inside the 
given intervals, the expressions tend to become unstable for scales outside 
of the chosen ranges. In table [| we have confronted the NLO and NNLO 
corrections to the LO. The NNLO corrections turn out not to be so large for 
values of [i so ft > 1.2 GeV. But this may happen accidentally since there are 
cancellations for the combination of scales used in this analysis. Therefore 
this should not be taken as an indication for small higher order corrections. In 
order to obtain a more conservative estimate, thus we have used the variation 
of the scales in our error analysis. 

A significant uncertainty also comes from Aq C d- By choosing Aq C d = 
330 ± 30 MeV we get an error of AM C = 50 MeV. The variation of the con- 
tinuum threshold sq and the error from the experimentally measured decay 
widths only have a small influence on the mass. We have summarised the 
results in table f|. Adding the errors in quadrature we obtain the charm pole 
mass: 

M c = 1.70 ± 0.13 GeV. (12) 
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Figure 3: Thick solid line: central PS-mass; thin solid lines: nips for 
Hsoft — 1-05 and 1.5 GeV; dashed lines: m PS for /ij ac = 1.2 and 1.7 GeV; 
dotted lines: mps for fihard = 1.4 and 2.5 GeV. 



Using the three-loop relation between the pole and the MS-mass which has 
been calculated recently []3~I], |3"^] , we obtain m c {m c ) = 1.20 ± 0.13 GeV for 
the MS-mass. In addition, we get an error from the perturbative relation 
between the two masses which does not converge for the charm case and can 
be of order Aqcd- By using the PS-mass, in the next section we will derive 
a more precise value for the MS-mass. 



6.3 Potential-subtracted mass scheme 

First, we have to choose a value for the separation scale fj, sep . This scale 
should be taken large enough in order to guarantee a perturbative relation 
to the MS-mass. On the other hand, it should be smaller than Mv. Both 
conditions cannot be well fulfilled at the same time. As a compromise value 
we will choose n sep = 1.0 ± 0.2 GeV. In this scheme, the pole contributions 
from the Greens function turn out to be smaller than in the pole mass scheme. 
The contributions from the condensates get more important here and we 
shall restrict our analysis to n < 6 to where these corrections are under good 
control. 

Using the same central values for the scales (|l(]) we obtain mps(^ S e P = 
1.0) = 1.35 GeV and from this value a MS-mass of m c {m c ) = 1.23 GeV. 
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n 


3 


4 


5 


6 


KA Poles 


0.59 


0.40 


0.27 


0.19 




1.16 


0.83 


0.63 


0.49 


KA N RResum 
J l n 


0.33 


0.32 


0.30 


0.29 


KA Counter 
J l n 


0.53 


0.47 


0.42 


0.37 


KA Continuum 
J l n 


0.42 


0.20 


0.09 


0.03 


KA Condensates 
J l n 


-0.07 


-0.09 


-0.12 


-0.14 



Table 5: Moments for different n with the parameters [i so ft — 1-2 GeV, 
fifac = 1.45 GeV, fi hard = 1.7 GeV, ^/s^ = 3.8 GeV and v sep = 0.3. 



f^sep 


0.8 


0.9 


1.0 


1.1 


1.2 


m PS 


1.38 


1.36 


1.35 


1.34 


1.33 


m c (m c ) 


1.20 


1.22 


1.23 


1.24 


1.26 



Table 6: Change of the masses for different values of id sep . 

The introduction of the intermediate mass definition leads to a reduced scale 
dependence: 

1.05 GeV < [i soft < 1.5 GeV : Am PS = 60 MeV , 

1.2 GeV < fi fac < 1.7 GeV : Am PS = 40 MeV , 

1.4 GeV < fx hard < 2.5 GeV : Am PS = 30 MeV . (13) 

In table [5] we again have collected the moments for different values of n. 
The influence of the pole contributions is clearly reduced and becomes even 
smaller for higher values of n. Since we have chosen the separation scale 
relatively large in comparison to Mv, the PS-mass lies a good deal lower 
than the pole mass and we are further away from the threshold region. The 
convergence of the perturbative expansion is improved in this scheme as well. 
In tabled we have varied the separation scale n sep between 0.8 GeV < [i sep < 
1.2 GeV. The error on the MS-mass is about 30 MeV. The change of the two 
masses is not directly correlated since the variation of fi sep also changes the 
relation between the masses. When compared to the pole mass scheme, the 
uncertainty from a s is reduced and the influence of the condensates becomes 
more important here. 

In table [7| we have listed the individual contributions to the error of mps 
and m c . Finally, we obtain the masses: 

m PS (fx sep = 1.0) = 1.35 ± 0.09 GeV , 

m c (m c ) = 1.23 ±0.09 GeV. (14) 
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Source 


Am PS 


Am c 


Variation of fJ, so ft 


60 MeV 


50 MeV 


Variation of /i/ ac 


40 MeV 


40 MeV 


Variation of Hhard 


30 MeV 


30 MeV 


Variation of fi sep 


30 MeV 


30 MeV 


Threshold s 


30 MeV 


30 MeV 


Experimental error 


10 MeV 


10 MeV 


Condensates 


20 MeV 


20 MeV 


Variation of v sep 


10 MeV 


10 MeV 


Variation of Aqcd 


10 MeV 


20 MeV 


Total error 


90 MeV 


90 MeV 



Table 7: Single contributions to the error of mps and m c . 



When we compare this value of the MS-mass to the one from the last section, 
Tn c (m c ) = 1.20 GeV, we see that the central value remains almost unchanged. 
This is not self-evident, as the dominating pole contributions are reduced in 
the PS-scheme and the relative influence of the individual contributions is 
shifted. 



7 Conclusions 

To summarise, our values obtained for the charm quark pole- and MS-mass 
follows: 

M c = 1.70 ±0.13 GeV, 
m c (m c ) = 1.23 ±0.09 GeV. (15) 

The obtained value for the pole mass lies somewhat higher than in previous 
sum rule analyses |28|,|29|]. In |28| the authors used perturbation theory to 



NLO resulting in a value of M c = 1.46±0.07 GeV. In the second investigation 
p9| , |33"H the analysis has been performed in the MS-scheme with perturbation 
theory to NLO. Using the NLO relation to the pole mass the author obtains 
m c {m c ) = 1.26 ± 0.05 GeV and M c = 1.42 ± 0.03 GeV. The author has 
also performed an analysis using resummation in LO with a value of M c = 
1.45 ±0.07 GeV. 

In our analysis the increased value of the pole mass is essentially due to 
large Coulomb contributions which have not been included in former analy- 
ses. As a consequence, the error becomes larger as well. In a recent analysis, 
the pole mass has been estimated from the charmonium ground state at 
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NNLO Here the authors obtained a pole mass of M c = 1.88±qH GeV. 
This value resulted from large corrections of the Coulomb potential to the 
ground state. 

During the last years, several lattice analyses obtained the following val- 
ues for the MS-mass: 



m c [m c ) 
m c (m c ) 
m c (m c ) 
m c (m c ) 



1.59 ±0.28 GeV [34] 

1.33 ±0.08 GeV ||] 

1.73 ±0.26 GeV || 

1.22 ±0.05 GeV J37], 



(16) 



Whereas the results from [35,37,38] are in good agreement with this anal- 
ysis, the investigations from [3^] and |36j obtain higher masses. For the 
time-being, the results are not conclusive and future lattice calculations for 
the charm mass might be of interest. Since other methods reveal significant 
uncertainties in the determination of the quark masses, QCD sum rules re- 
main amongst the most precise tools to extract these fundamental quantities. 
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